1 



The Clusters AgeS Experiment (CASE): 
The blue straggler Star M55-V60 r^aught amidst 
rapid mass exchangebl 



M. R o z y c z k a^, J. K a 1 u z n y^, I. B. T h o m p s o n^, 
S. M. R u c i n s k i^, W. P y c and W. Krzeminski^ 

-'^Nicolaus Copernicus Astronomical Center, ul. Bartycka 18, 00-716 Warsaw, 

Poland 

e-mail: (jka, mnr, wk, "wp)@camk. edu.pl 
^The Observatories of the Carnegie Institution of Washington, 813 Santa 
Barbara Street, Pasadena, CA 91101, USA 
e-mail : ian@obs . carnegiescience . edu 
■^Department of Astronomy and Astrophysics, University of Toronto, 50 St. 
George Street, Toronto, ON M5S 3H4, Canada 
e-mail: rucinski@astro . utoronto . ca 



ABSTRACT 

We analyze light and velocity curves of the eclipsing blue straggler V60 in the field of the 
globular cluster M55. We derive A/p=1.259±0.025MQ , _Rp=1.102±0.021Rq , M^°'=3.03±0.09 
mag for the primary and Ms = 0.327 ± 0.017 , = 1.480 ± 0.011 /?0 , M^*""' = 4.18 ± 0.12 
mag for the secondary. We measure an apparent distance modulus {m — M)v = 14.04 ±0.09 
mag. Based on the systemic velocity, distance, and proper motion of V60 we conclude that 
the system is a member of the cluster and argue that its present state is a result of rapid but 
conservative mass exchange which the binary is still undergoing. We report a peculiar blue 
excess on the ascending branch of the primary eclipse of V60 and discuss its possible origin. 
globular clusters: individual: M55 - blue stragglers - stars: individual: V60-M55 



1 Introduction 

The eclipsing variable M55-V60 (henceforth V60) was discovered in 1997 within 
the CASE project (Kaluzny et al. 2010). Located 27 arcsec from the center 
of M55, it is a proper motion member of the cluster (Zloczewski et al. 2011). 
At Vmax = 16.8 mag and {B — V)max = 0.41 mag the system belongs to the 
population of blue stragglers. Kaluzny et al. (2010) compiled photometric 
data from several seasons and found that the light curve was characteristic of 
short-period low-mass Algols: deep primary eclipses with DV 1.8 mag were 
followed by shallow secondary ones with DV^i 0.20 mag. The orbital period, 
equal to 1.183 d, was found to systematically increase. Such properties indicate 
a semidetached binary in the mass-transfer phase, with the donor being the less 
massive (secondary) component filling its Roche lobe. 

In this paper we analyze the velocity curve of V60 and refine the preliminary 
photometric solution obtained by Kaluzny et al. (2010). Spectroscopic and 
photometric data are described in Sect. [2] and Sect. [3l The analysis of the data 
is detailed in Sect. |4] and our results are discussed in Sect. [S] 



*This paper uses data obtained with the Magellan 6.5 m telescopes located at Las Cam- 
panas Observatory, Chile. 
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2 Photometric observations 

The photometric data collected between 1997 and 2009 were discussed in Kaluzny 
et al. (2010), but for completeness we repeat here the most important points of 
their discussion. 

According to their findings, season-to-season variations of the light curve in 
quadratures and in the secondary eclipse do not exceed 0.010-0.015 mag in V- 
band. A significantly stronger variability is observed in the primary minimum, 
whose depth ranged from V = 18.62 mag in 1999 to V — 18.57 mag in 2007 
and 2009. The primary minimum is symmetric in V, allowing for a precise 
determination of times of minimum (however, as we show and discuss in Sect, 
m the symmetry is broken in B). The 0~C diagram of the times of minimum 
(see Fig 3. in Kaluzny et al. 2010) indicates that the orbital period of V60, which 
in 2008 was equal to 1.1830214±0.0000007 d, is increasing at a rate dP/dt = 
3.0 X 10^^; i.e. it would be doubled in just ^10^ years. 

A preliminary photometric solution favors a semi-detached configuration 
with the secondary filling its Roche lobe, consistent with the observed be- 
havior of the period. The luminosity ratios of the components in V and B 
bands (obtained from the same solution at quadratures) yield apparent colors 
{B — V)p ~ 0.3 mag and [B — V)s ~ 0.9 mag, respectively for the primary and 
the secondary. With Vp = 17.11 mag the primary of V60 is located among the 
SX Phoenicis pulsating variables on the color-magnitude diagram (CMD) of 
M55. A sinusoidal variation with an amplitude of ^0.005 mag and a period of 
0.03087±0. 00007 d is indeed observed in out-of eclipse sections of the nightly 
F-light curves. 

For the present analysis we use observations from five seasons in which both 
eclipses were covered (i.e. 1999, 2006, 2007, 2008 and 2009). The data were 
phased according to the times of minima given in Table 4 of Kaluzny et al. 
(2010) and merged into the composite B and V light curves (with 663 and 2666 
points, respectively) shown in Fig. [B 

3 Spectroscopic observations and orbital param- 
eters of V60 

Our radial velocity data are based on eight observations obtained with the MIKE 
Echelle spectrograph (Bernstein et al. 2003) on the Magellan II (Clay) telescope 
at the Las Canipanas Observatory. Seven of them were made in 2004 between 
June 27th and October 3rd, and the eight one in 2005 on September 11th. Each 
observation consisted of two 1200-1800 s exposures interlaced with an exposure 
of a Th/Ar lamp. For all observations a 0.7 x 5.0 arcsec slit was used, and 2x2 
pixel binning was applied. At 438 nm the resolution was ^2.7 pixels at a scale 
of 0.0043 nm/pixel. The spectra were processed using a pipeline developed by 
Dan Kelson following the formalism of Kelson (2003). In the blue channel of 
MIKE (380 - 500 nm) the average S/N ratio ranged between 21 and 25. As 
in the red channel (490 - 1000 nm) it was markedly worse, for radial velocity 
measurements data from the blue channel only were used. 

The original data had a high resolving power (about 40,000), but because 
the S/N per pixel was rather low, we rebinned the continuum- rectified spectra 
to a log(A) scale, convolved them with a Gaussian with FWIIM=15 km/s, and 
rebinned them to 3 km/s steps (i.e. 5 pixels per a resolution interval or 2.5 
times over-sampling) . This procedure improved the quality of the spectra and 
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adjusted the resolution to better reflect the actual rotational broadening of the 
lines. The final spectra have 22,300 data points in the wavelength interval of 
400 - 500 nm, where the S/N ratio was the highest. The data reprocessed in this 
way were analyzed with the help of a code based on the broadening function 
(BF) formalism of Rucinski (2002). The BFs were determined over 261 steps 
of 3 km/s in a process which is effectively a least squares solution of 22,300 
linear equations with a 22,300/261 — 85 times over-determinacy. Templates 
with [Fe/H] = — 2.0, almost exactly matching the value of -1.94 given by Harris 
(1996; 2010 edition), were selected from the Synthetic Stellar Library of Coelho 
el al. (2005). Measurement errors, estimated based on the residuals from the fit 
described in the next paragraph, amount to 3.6 km s~^ (primary) and 1.9 km s""'^ 
(secondary). Rotational profile fits which are performed automatically within 
the BF formalism yielded Vrotsini equal to 46.1 ±1.4 km s~^ for the primary 
and 63.2 ±2.2 km s~^ for the secondary (corrections for the assumed Gaussian 
smoothing with FWHM — 15 km/s were taken into account). 

The velocity curve was fitted with the help of the spectroscopic data solver 
written and kindly provided by Guillermo Torres. On the input to the solver 
the observed velocities had to be transformed into mass-center velocities of the 
components in order to account for their asphericity. This was done iteratively 
by applying the solver to the uncorrected velocities, feeding the solution into 
PH0EBE31a implementation (Prsa & Zwitter 2005) of the Wilson-Devinney 
model (Wilson & Devinney 1971; Wilson 1979), finding the corrections, and 
applying the solver again to the corrected velocities. The advantage of this 
procedure is that Torres's code automatically calculates the errors of the fitted 
orbital parameters, which otherwise would have to be estimated with PHOEBE 
by Monte-Carlo techniques. The measured orbital velocities are listed in Ta- 
ble [1] together with mass-center corrections and residuals from the fit. The 
fitted velocity curve generated by PHOEBE is shown in Fig. [2l and the orbital 
parameters obtained from the fit are listed in Table [2] together with formal 1-cr 
errors returned by the fitting routine. 

4 Photometric solution and system parameters 

A closer look at the primary eclipse in B reveals a significant asymmetry iden- 
tifiable in each observing season: the ascending branch is by up to 0.2 mag 
bluer than the descending one (see Fig. [3]). As such an excess is not possible 
to model within the standard approach, we were forced to discard the affected 
observational points between phases 0.0 and 0.07. In principle, the pulsational 
modulation with an amplitude of~0.005 mag in V mentioned in Sect. [2l which 
is at least partly responsible for the vertical scatter of points in Fig. [1] should 
also be removed. However, an amplitude this low is suggestive of nonradial pul- 
sations which are commonly observed in SX Phe stars (e.g. Pych et al. 2001; 
Olech et al. 2005), and the results of such a procedure applied to the primary 
eclipse would be unreliable. 

The photometric solution was also found with the help of PHOEBE utility 
which employs the Roche geometry to approximate the shapes of the stars, 
uses Kurucz model atmospheres, treats reflection effects in detail, and, most 
importantly, allows for the simultaneous analysis of B and V data. For the 
radiative envelope of the primary we adopted a gravity brightening coefficient 
gp=1.0 and a bolometric albedo Ap— 1.0. The same coefficients for the convective 
envelope of the secondary were set to gs — 0.32 and Ag — 0.5. The effects of 
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reflection were included. Limb darkening coefficients were interpolated from the 
tables of Claret (2000) with the help of the JKTLD codeQ Full synchronization 
of both components was assumed. 

An approximate temperature of the primary Tp was calculated from the 
dereddened B — V index obtained by Kaluzny et al. (2010), using a color- 
temperature calibration based on the data from the Dartmouth Stellar Evo- 
lution Database (Dotter et al. 2008). We decided to employ the synthetic 
calibration because the starting value {B — V)po = 0.23 mag was too close to 
the applicability limit of the empirical calibration compiled by Casagrande et 
al. (2010) which is formally valid for 0.18 mag <{B — V)< 1.29 mag (this issue 
is further discussed in the last paragraph of this Section). For dereddening a 
value -^(i? — V^) =0.08 mag with an assumed error of 0.01 mag was used (Harris 
1996, 2010 edition). The next approximation was found based on PHOEBE- 
provided contributions of each component to the total light at quadratures in 
B and V bands which allow to calculate the updated observed {B — V)p. The 
updated [B — V)p was dereddened and translated into temperature the same 
way as before. The procedure was repeated until the convergence was reached 
at Tp = 8160 ±140 K (the error is due to uncertainties in calibration, reddening 
and zero points of B and V photometry). The temperature of the secondary, 
Tg, was automatically adjusted by PHOEBE. 

The fitted photometric parameters of V60 are listed in Table [3] together with 
the errors estimated based on additional fits to each of the seasonal B and V light 
curves (the ascending branch of the primary eclipse was removed from each in- 
curve while fitting) . Table |4] contains the final absolute parameters of the system, 
hereafter referred to as the standard solution. The residuals from the final fit 
are shown in Fig. IH The blue excess on the ascending branch of the main eclipse 
is clearly visible; it is also evident that the largest T^-residuals occur within the 
main eclipse. Apparently, the stream(s) of gas between the components generate 
an additional light in the system which PHOEBE is not able to account for. The 
likely variability of such a light source might be responsible for the enhanced 
scatter observed in y-band at the bottom of the main eclipse. 

For rp=8160 K and the primary's gravitational acceleration gp—AA obtained 
from our solution, the bolometric correction amounts to -0.07 mag (Dotter et al. 
2008; on-line version of Dartmouth Stellar Evolution Database). With Mfco/ = 
3.03 ±0.09 mag (Table 4), the absolute V magnitude of the primary is Mvp = 
3.10±0.09 mag. From the light curve solution we obtained Vp — 17.144±0.014 
mag, where the error includes uncertainty of the zero point of the photometry. 
The apparent distance modulus of V60 is then (m— M)y = 14.04±0.09 mag. This 
is consistent with the value of 13.89 listed by Harris (1996; 2010 edition). With 
E{B — V) = 0.08±0.01 mag, the foreground absorption in T^-band is Ay = 0.25± 
0.03 mag, corresponding to an absolute distance modulus (m— M)o = 13.79±0.10 
mag. A very similar value follows from the result obtained by isochrone fitting 
(Dotter et. al. 2010), who derived an apparent distance modulus of 13.88 
mag for the ACS/F814W filter. Since from the calibration of Girardi et al. 
(2008) one gets Apsi4w — 0A5 mag, the absolute distance modulus becomes 
(m-M)o = 13.73 mag. 

We would like to note, however, that the primary's bolometric luminosity 
we derive strongly depends on the effective temperature deduced from the color 
index. If instead of rp=8160 K we used rp = 7450 K resulting from the empirical 
calibration of Casagrande et al. (2010), we would obtain an unacceptably low 



^Written by John Southworth and available at www.astro.keele.ac.uk/jkt/codes/jktld.html 
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value of the absolute distance modulus (m — M) = 13.39 mag. This casts come 
doubts on the validity of that calibration for blue low-metallicity stars, and 
indeed, a closer inspection of Fig. 14 in Casagrande et al. (2010) strongly 
suggests that for [Fe/H] = -2.0 and B-V < 0.35 mag their T^ff -~{B-V) 
relation is but an extrapolation. 

5 Discussion 

The important conclusion following from our analysis is that the mass transfer 
in V60 must be to a good approximation conservative. This is because the total 
mass of the system, M = 1.59 Mq, is almost exactly two times larger than the 
M55 turnoff mass of ■^O.S M0 (e.g. Zaggia et al. 1997). This, together with 
the rate of period lengthening and system parameters found in Sect. |4l implies 
an orbital expansion rate dA/dt = 3.4:X 10~^ Roy~^ and a mass transfer rate 
dM/dt = 1.4 X 10~^ M0y~^. Thus, V60 is in a phase of a rapid mass exchange, 
which, given the low mass of the H-shell burning secondary, cannot last longer 
than a few hundred thousand years. With this in mind, it is interesting to see 
that the present primary must be quite far from thermal equilibrium, as its 
temperature is over a thousand K lower than the temperature of a 1.1 Rq star 
on the 1.26 Mq evolutionary track, which according to Dotter et al. (2008) 
should be equal to 9350 K. We also note that the original mass ratio must have 
been close to unity - otherwise the original primary would have left the main 
sequence much earlier, and the system would be in a more advanced evolutionary 
state (e.g. detached, or semidetached with the present primary filling its Roche 
lobe). 

The blue excess of the system between phases 0.0 and 0.07 (hereafter BEx) 
is a puzzling feature which to our knowledge has no counterpart in other Algols. 
Suspecting it might be spurious, we checked if there was a systematic shift 
between phases calculated for B and V data points. None was detected. We 
then tried various means to remove BEx or at least to make it weaker. First, we 
fitted the complete B curve allowing PHOEBE to find a global phase shift A(p 
that would center the main B-minimum at phase 0. Subsequently, with Aip kept 
fixed, both B and V curves were fitted. This procedure indeed made the blue 
excess smaller, but the V-fit became significantly worse both in the primary and 
the secondary minimum. In a second series of experiments we fitted the curves 
with the descending branch removed from either B or B and V primary minima, 
allowing for a global phase shift or keeping it fixed at 0. Since the residuals did 
not look any better in any of those fits, we had to accept the standard solution 
as the best one we were able to obtain. The encouraging finding was that the 
parameters of all the trial fits did not differ from the standard ones by more 
than a factor of ~1.5 times the errors given in TableU) In particular, this means 
that the masses of the components are determined sufficiently accurately for 
the conclusion concerning the conservative mass transfer to be firmly footed. 
We note, however, that while the mass exchange was conservative or nearly 
conservative, the system must have lost an appreciable amount of the orbital 
angular momentum Jorb, possibly by magnetic braking (see e.g. Eggleton and 
Kiseleva-Eggleton 2002). At the present value of Jorb the minimum orbital 
separation (which is achieved when Ms = Mp) would be 0.94 R© - far too small 
to accommodate two 0.8 Mq stars. Of course, the loss of Jorb must have occurred 
after the equalization of masses, i.e. when the original mass ratio was reversed 
and the orbital separation began to increase. 
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Numerous spectroscopic and photometric effects observed in Algol-like sys- 
tems are attributed to three dynamical agents: a stream of matter flowing out of 
the secondary, a regular or transient accretion disk formed by the stream around 
the primary, and shock waves excited where the stream hits the disk or directly 
impacts the primary. V60 is a compact system in which there is no space for 
a regular accretion disc (e.g. Richards 1992), so that the stream directly hits 
the primary and induces shock waves(s) in its atmosphere. BEx, if real, must 
also originate in a shock. The problem is that the shock cannot be localized in 
the standard place where the trailing hemisphere of the secondary is hit by the 
stream (SI in Fig. [5]): if that were the case, then BEx would show up on the 
descending branch of the main minimum. We calculate that if a particle leaving 
the Li point is to miss the trailing hemisphere and land onto the leading one 
(path a ending at S2 in Fig. [5]) an initial velocity of ~0.5Vorb is needed, where 
Vorb is the orbital velocity of the secondary around the mass center. Since in 
reality the initial velocity is on the order of the thermal velocity, i.e. about a 
tenth of Vorb, the stream must follow path b in Fig. [S] and hit the star at SI. 
Another shock at 82 can only be formed by the matter reflected off the primary 
and flowing along the continuous line in Fig. [5] Such an interpretation is not a 
new one - in fact, it closely resembles the schematic Algol model derived from 
spectroscopic observations by Gillet et al. (1989) and shown in their Fig. 3. 
Unfortunately, V60 is too faint for a detailed spectroscopic study, and the only 
observational proof for the existence of S2 would be a UV emission coinciding 
in phase with BEx. 

We note that speculations about the possibility of stream reflection are not 
entirely unrealistic: a partial reflection was observed in models of flows in a 
cataclysmic binary obtained by Rozyczka (1987). Admittedly, in those simula- 
tions the stream reflected off the accretion disk, but the physics involved was 
sufficiently simple that the results may be safely applied to a stellar atmosphere. 
It should also be mentioned that a reflection from the primary component in an 
Algol-type system was observed in simulations performed by Richards (1998), 
however it was to a large extent predetermined by boundary conditions. 

Assuming the above explanation of the origin of BEx is feasible we still have 
to explain why the emission from S2 is only visible in a narrow range of phases, 
and why the emission from SI is not visible at all. The only possibility we see 
is that both shocks are hidden behind gas stream(s): SI entirely, and S2 partly. 
This speculation can be verified by detailed hydrodynamical simulations which 
are beyond the scope of this paper. 
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Table 1: Radial velocity observations of V60. 



HJD 


Phase 


Vp 




dvp 


dvs 


{0-C)p 


{0-C)s 


-2453000 




[km/s] 


[km/s] 


[km/s] 


[km/s] 


[km/s] 


[km/s] 


183.79260 


0.7523 


221.66 


-8.87 


0.07 


-3.49 


0.27 


0.44 


206.73462 


0.1451 


133.18 


319.61 


-0.07 


1.72 


-1.84 


1.21 


272.53145 


0.7626 


226.01 


-7.92 


0.07 


-3.46 


4.77 


0.85 


275.56501 


0.3269 


129.41 


337.27 


-0.03 


1.23 


-0.98 


0.70 


276.57721 


0.1825 


126.29 


336.59 


-0.07 


2.54 


-2.88 


-3.51 


281.51625 


0.3574 


138.33 


317.70 


-0.02 


-0.66 


2.89 


-1.28 


282.51218 


0.1994 


134.86 


348.79 


-0.07 


2.92 


7.54 


1.94 


989.57690 


0.8788 


206.18 


41.17 


0.06 


-1.34 


-0.78 


-2.95 



Columns dVp and dVs contain corrections which have to be added to velocities 
Vp and Vs in order to transform them into mass-center velocities of the primary 
{p) and the secondary (s). 



Table 2: Orbital parameters of V60 



Parameter Unit Value 

~P d 1.1830214±0.0000007'' 

7 km s-i 173.14±0.68 

Kp km s-i 48.32±1.59 

K, km s-i 185.96±0.87 

e O.O'' 

(Tp km 3.56 

fjs km s~^ 1.91 
Derived quantities: 

^sini Rq 5.476±0.046 

MpSin^z Mq 1.251±0.025 

Mssin^i Mq 0.325±0.017 

^Obtained from photometry 
Assumed in fit 
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Table 3: Photometric parameters of V60. 



Parameter 


Unit 


Value 


i 


deg 


86.4±0.5 


e 




Qa 


T 

-i-p 


K 


8160'" 


T 


K 


5400±77 


Rp 


Rq 


1.102±0.019 


Rs 


Rq 


1.480" 


{Lp)v 


% 


76.5±1.8 


{Lp)B 


% 


85.3±1.5 


(JV 


mmag 


18 


CTB 


mmag 


12 


'^Assumed in 


fit 





Obtained from the synthetic {B — V)—T calibration based on the Dartmouth 
Stellar Evolution Database 
^For fixed A = 5.487 Rq 



Table 4: Absolute parameters of V60 



Parameter 


Unit 


Value 


P 


d 


1.1830214± 0.0000007 


A 


Rq 


5.487±0.047 


i 


deg 


86.4±0.5 


e 




ga 


Mp 


Mq 


1.259±0.025 


M, 


Mq 


0.327±0.017 


Rp 


Rq 


1.102±0.021'" 


Rs 


Rq 


1.480±0.01l'' 


T 


K 


8160±140'= 


Ts 


K 


5400±160'* 




mag 


3.03±0.09 


jjbol 


mag 


4.18±0.12 


"^Assumed in 


fit 





^Errors include uncertainty in A 

'^Combined error of calibration, reddening and zero points of B and V photom- 
etry 

'^Error includes uncertainty in Tp 
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Figure 1: Light curve of V60 in bands V (upper) and B (lower). Data from five 
observing seasons are shown. 
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Figure 2: Radial velocity curve of V60 (up er panel), and residuals from the fit 
(lower panel). 



12 




Figure 3; Asymmetry of the primary eclipse of V60 in the _B-band. Points: 
observational data from five observing seasons phased as explained in Sect. [2j 
Crosses: descending branch reflected in phase. 
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Figure 4: Residuals from the fits to V and B light curves. The i3-residuals in- 
clude points on the ascending branch of the primary eclipse which were omitted 
during the fitting (see text for explanations). The vertical scale is the same in 
both panels. 
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Figure 5: Schematic view of V60 with sizes of the stars and binary separation 
drawn to scale. See text for explanations. 



